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A treatment of edge magnetoplasmons (EMP), based on a microscopic evaluation of the local 
contributions to the current density, is presented. It is valid in the quantum Hall regime for filling 
factor = 1 or 2 and low temperatures when the dissipation is localized near the edge. The confin- 
ing potential, flat in the interior of the channel, is assumed smooth on the magnetic length scale 
but sufficiently steep at the edges that the density profile is sharp and the dissipation considered re- 
sults only from electron intraedge-intralevel transitions due to scattering by piezoelectrical phonons. 
For wide channels there exist independent EMP modes spatially symmetric and antisymmetric with 
respect to the edge. Certain of these modes can propagate nearly undamoped even when the dissipa- 
tion is strong and are thus termed edge helicons. In contrast with well-known results for a spatially 
homogeneous dissipation within the channel, we obtain that the damping of the fundamental EMP 
^ ■ is not quantized and varies as T'^ or T^'^, where T is the temperature, in the high- and low-frequency 

' limits, respectively. The characteristic length of the resulting dispersion relation and of the charge 

density distortion is Iq. The screening of the metallic gates, when present, is taken into account. 
PACS 73.20.Dx, 73.40.Hm 
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I. INTRODUCTION. 



In the past few years there has been considerable interest in edge magnetoplasmons (EMP) as well as in other edge 
. excitations of two-dimensional (2D) electron systems in the presence of a magnetic field B |^]- For a 2D system 
^ ■ with a vertical conductivity drop at the boundaries, it has been shown that the dissipation can determine the EMP 
O _ dispersion relation and the spatial structure in an essential manner even in the regime of the quantum Hall effect 
I ^ I (QHE). In this work the properties of the EMP are expressed in terms of the components of the magnetoconductivity 
tensor of an infinite 2D system. Moreover, due the very low frequency uj of the EMP. the dispersion relation could be 
written in terms of the static magnetoconductivity tensor. 

The distance of the "center of gravity" of the EMP charge from the edge, which coincides with the characteristic 
length over which the transverse to the edge electric field Ey of the EMP decreases, is given, for \kxiv\ Ij by [|] 



> 

oo 

oo . 

On 1 Here ayy and ayx are the conductivity components of an infinite 2D system, is the EMP wave vector, and iy denotes 
a characteristic length determined by Eq. (10) of Ref. 0]. In the QHE regime for typically observed EMP, Eq. (|l|) 
^ , gives ivc ^ 1/im. Now, in Ref. [ p^ we have shown theoretically, and in agreement with experimental observations, 
that in the QHE regime for sufficiently smooth confinement the dissipation is due to intralevel-intraedge transitions 
of electrons scattered by piezoelectrical phonons and occurs mainly near the edges of channels. In the linear response 
regime this is the main dissipation for channels of width W ^ 100 /im and temperatures T ^ IK if the group velocity 
Q of edge states, Vg, is larger than the speed of sound s. As for dissipation in the bulk, it is exponentially suppressed 
Q . for T ^ 0. Given that and the fact that the dissipation, when deriving Eq. (|l|), occurs in the bulk, we expect the 
L| ' properties of the EMP to be strongly modified when the dissipation is locahzed near the edges. 

^ i The above expectation is further supported by the results of Ref. which pertain to EMP for a smooth, un- 
■ perturbed electron density profile which contrasts sharply with that of Ref. where the density drops vertically 

" ~ ~ '^1 acoustic EMPs were obtained in Ref. |l^. Further, our results 
. Ql and |ll|] , show significant differences from them even when the 
dissipation is very weak. For the very low temperatures that we consider here, fc^T <g; hvg/io, and the assumed 
smooth confining potential on the scale of (q {vg > s), the unperturbed electron density no{y), normalized to the bulk 
value no, drops essentially, on the scale of £o, only near the edge. More explicitly, for the potential that is specified 
at the beginning of Sec. II, we calculate no{y)/no = {1 -|- ^[{yre — y)/^o]}/2, where yre is the coordinate of the right 
edge and <&(?/) the probability integral. In Fig. 1 we show this calculated density profile (short-dashed curve) together 
with those assumed in Refs. (solid curve) and |ll| (long-dashed curve). The profile of Ref. is obtained with 
no(y)/rio = {2/'K)arctany^ {yre ~y)/a and a/io = 20; it corresponds approximately to a = 2000A. As can be seen, 



5h ] at the edges. In addition to the modes of Ref. 
- - 1 despite their partial similarity with those of Refs. 
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the three density profiles are very different from each other. As will be shown in this paper, combining our density 
profile with the localization of the dissipation near the edge leads to strong medications of the EMPs results. 

These modifications, as well as new EMP resulting from the microscopic treatment of the problem, are the subject 
of this work. The description of the inhomogeneous current density in the quasi-static regime is carried out using 
the results of Ref. jl^. We consider the QHE regime, mainly u — 1 and partly v — 2, for samples with sufhciently 
large in- plane dimensions, as is typical in EMP experiments, that inter-edge electron transitions and the inter-edge 
Coulomb interaction can be neglected. 

In Sec. II we start with the expressions for the inhomogeneous current densities and conductivities and derive the 
integral equation for EMP. In Sec. Ill we derive the dispersion relation for very low temperatures and in Sec. IV 
we describe in detail the new edge waves. Finally, in Sec. V we compare our theory with the experiment and make 
concluding remarks. 

II. BASIC RELATIONS 

A. Inhomogeneous current density in quasi-static regime 

We consider a two-dimensional electron gas (2DEG), of width W, of length — L, and of thickness zero, in the 
presence of a strong magnetic field B parallel to the z axis. The 2DEG is confined along the y axis. For simplicity 
we take the confining potential as parabolic at the edges: Vy = 0, for yi < y < y,-, Vy = m*n'^{y — yrY /2 for 

y > yr > 0, and Vy = m*^'^{y — yiY /2 for y < yi < 0. Because in real EMP experiments W > 0.1cm, we can 
assume [W — {yr — yi)]/W <C 1. Moreover, we will assume that ^ 1 such that it is possible to consider 

an EMP along the right edge of the channel, of the form A(oj,kx,y)exp[—i{ujt — k^x)], totally independent of the 
left edge. We consider only linear responses. For definiteness, we take the background dielectric constant e to be 
spatially homogeneous. We consider B strong enough that only the n — Landau level (LL) is occupied. For the 
ly = 2 QHE regime we will neglect the spin-splitting. As for the v — 1 QHE, we will assume that the spin-splitting, 
caused by many-body effects, is strong enough to neglect the contribution related with the upper spin-split LL. We 
assume a lateral confinement smooth on the scale of the magnetic length Iq — (?i/m*a;c)^^^ such that Vl lUc, where 
ujc = \e\B/m* is the cyclotron frequency. Further, we will approximate u! — (u)^ + Sl^)^/^ by lUc- 

Because the EMP is practically quasi-static and its wavelength A ^ 1cm is very large, we expect, in analogy with 
well-known results that follow from Maxwell's equations ||l^, the associated electric field Ex{x, y, t) to have a smooth 
dependence on y on the scale of io = {h/m*u!cy^^ , i. e., Ex{x, y, t) = E^iy) exp[— — kxx)]. Physically it is clear: 
the dependence of Ex{x,y,t) on y, as expressed through Maxwell's equations, is related to that on x which has a 
characteristic scale A. Thus, E^iy) should have the same scale A and be definitely smooth on the £q scale. This is a 
general result and applies to the case treated in Ref. j|] . Then using the results of Ref. iQ , we obtain the components 
of the current density in the form 

Jy iy) = '^yy {y)Ey (y) + a^], {y)Ex (y), (2) 

jx{y) = axx{y)Ex{y) - a°x{y)Ey{y) + vgp{uj,kx,y). (3) 

Here we have suppressed the exponential factor exjp[—i{ujt — k^x)] common to all terms in Eqs. (^) and (^. It is 
understood that E^{y) depends on cu and k^. As follows from Refs. and [Q, <Jyy{y) is strongly (exponentially) 
localized at the edge, within a distance ^ £o from it, for hvg ^ ksTlo. The last term on the riht-hand side (RHS) 
of Eq. (^), absent in Ref. [|4j, represents a convection contribution to the current density along x, associated with 
the wave, and is due to a distortion Sp of the charge localized near the edge; we denote it by p(Lu,kx,y) in order 
to simplicity the notation as it occurs frequently. Notice that in Ref. the contributions to the components of 
the current density are microscopically obtained for the electric field components smooth on the scale of Iq. This 
condition holds for the contributions oc E^iy) in Eqs. (||) and (^) but is not well justified for those cx Ey{y). We 
assume that the latter can be reasonably approximated by those obtained microscopically when Ey(y) is smooth on 
the scale of io- The assumption is equivalent to neglecting possible nonlocal contributions to the current density 
oc / dy'afiy{y,y')Ey{y'); then it follows that axy'^ = —ayx'^. For = 1 we have jlj] 

'^lx{y) = 7rT dyoafao'i'liy-yoa), (4) 

where a = {0,kxa}, yoa = ^okxa, ^n(y) is a harmonic oscillator function, and fao = fo{kxa) = + exp{{Eao — 
Epo) /ksT)] is the Fermi-Dirac function. Epo is the Fermi level counted from the bottom of the lowest electric 
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subband. For T = and near the right edge we have cr°^(j/) = (e^/47r?i){l + ^[{yre — y)/^o]}, where $(a;) is the 
probabihty integral, yre — elkre, and /o(fcre) = 1/2. That is, a°{y) near the edge decreases on the scale of £o and 
behaves as the density of the short-dashed curve of Fig. 1 . Considering only the right edge and the flat part of the 
confining potential, for yi < yoa < yr we have Eao = hLUc/2 and for j/qq > yr we obtain 

E^o = Eo{k,^) = hLu,/2 + m*n^{yoc. - yrf/2. (5) 

We consider only the interaction of electrons with phonons, and neglect that with impurities, since the former is the 
most essential for the assumed conditions Because of the very smooth dependence of Ex{y) on £o, we can assume 
that (Jxxiy) can be approximated by (Jyyiy) which follows from Eq. (16) of Ref. |jl^ as 

^yy^'^^ " AhLkBT ^ |Cqp9x[/o(fcxQ - qx) - /o(fcxa)]<5[£^o(fcxQ) - Eo^k^a - Qx) - ^q] 

fcxaq 

xe-(«^+«2)^o/2 sinh-^^) [^l{y - yoik,^ - Qx)) + ^Uv - yo{kxo.))]. (6) 
For the low temperatures pertinent to the quantum Hall effect we consider only the standard acoustical (DA-) or 



piezoelectrical (PA-) phonons for which ujq = sq, where s is the speed of sound, and q = ylx Qy ^ Qz- Then 
|Cqp = (c / LxLyLz)q^^ , where +1 is for DA- and —1 for PA-phonons, respectively. 



B. Integral equation for EMP with dissipation at the edges 

Using Eqs. (^-(^), and (^), we can write the continuity equation, linearized in 6p{uj, kx,y) = p{i~o, kx, y), as 
- i{uj - kxVg)p{uj, kx,y) + ik (y) Ex{uj,kx,y) - (Jyx{y)Ey{uj,kx,y)] 



+ -g:^[<^yy{y)Ey{uj,kx,y) + a°^{y)Ex{LJ,kx,y)] =0. (7) 



In terms of the potential 4){u!,kx,y) the electric field components are Ex{uJ,kx,y) = —ikx4'{w,kx,y) and 
Ey{uj,kx,y) = --§-<t){uj,kx,y). Then Eq. (|) gives 



- i{uj - kxVg)p{uj, kx, y)+klaxxiy)4'i^, kx, y) 
Now using Poisson's equation we obtain 

/•OO 



Wyy{y)jr(l^{^^kx,y)] - ikx(j){uj,kx,y)—(jlx{y) = 0- (8) 



^,kx,y)^-l dy'Kn{\kx\\y~y'\)p{uj,kx,y'), (9) 



where Kq{x) is the mofified Bessel function; 4> and p pertain to the 2D plane. From Eqs. and we obtain the 
following integral equation for p{lu, kx, y) 

- i{uj - kxVg)p{uj,kx,y) + ^{klaxxiy) - ^kx-^Wlxiy)] 

-^vv^y^^ - di^Mv)]^} J dy'Ko{\kx\\y~y'\)p{uj,kx,y')^0. (10) 

The value of (Tyy{y) is significantly different than zero only near the edges of the channel. The same holds for the 
values of (Jxx{y) and of dayx{y) / dy; for Hvg '3> ioksT this can be seen from Eqs. (^ and (^ which show that 
(Jyyijj) and doyxiy) / dy are exponentially localized within a distance ~ from the right edge at yre = yr + ^yr- 
We have Ayr = igke where fcg = {ujc/1iil.)y/2m*Ap is the characteristic wave vector associated with an edge state, 
Ap = Epo — fii^c/2, and W — 2yre- For kxa = kre = yrl^ + we have fo{kre) = 1/2 and 



_ I dEoikre) _ hn^ke _ 2Af Vt 
^ h dk^n m*ujl: V m* ojr' 
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Eq. (HH) can also be written as Vg = Ee/B, where i?e = VL\/2m* /S.p /\e\ is the electric field describing the influence 
of the confining potential. For a dissipationless, 2D classical electron liquid we have, for finite w, (Jyy(ij) = axxjy) — 
ie'^nQ{y)uj /m* {lo^ — uj"^) and cr°^(y) = —e'^no{y)uJc/m*{uj'^ — w^), where hq is the electron density. Then Eq. {Wij 
becomes identical with Eq. (4) of Ref . . In addition, if we assume that the conductivity components in Eq. ( ^) 
are independent of y, for |y| < W/2, and <Jyy{y) — (Jxxiy), Eq. (|lO|) takes the form of Eq. (15) of Ref. after 
integration over z. 

Equations and ( |lO|) apply to a 2DEG in the absence of metallic gates. Sometimes a metallic gate is placed on 
the top of the sampl e [tLOj at a distance d from the 2DEG. As shown in the Appendix, for a gated sample the kernel 
i^o in Eqs. and (|lO|) .s replaced by Rg — A'o(|fcj,||y — y'\) — K(){\kx\\/ {y — y'Y + AdP). If this gate is replaced by 
air, then is replaced by i?, - i^o(|fc.||y - y'\) + [(e - l)/(e + l)]K^{\kxW {y ^ y' f + A(P). 



III. EMP DISPERSION RELATION 

We consider very low temperatures that satisfy the inequality TiVg ^ ioksT. From Eqs. and it follows 
that day^{y)/dy = — (e^/27r?i)^'o(?/ — ?/re); also, g-yy jy) and <Jxx{y) behave as 5'o(y — J/re) and hence are strongly 
concentrated near the edge. It follows from Eq. ( |lO| ) that p{ijj,kx,y) is also concentrated near the edge. Integrating 
Eq. ( ]lO| ) over y, from y^e ~ to y^e + ^2/ with Ay ^ to, we obtain 

/oc 
dy'^liy' -yre)Ko{\kx\\y~v'\)]^0; (12) 
-oo 

here S* = (2/e){—ik'^axx + ^^0"^^), cr^^ — e^/2Trh is the Hall conductivity in the bulk as follows from Eq. (|^), and 
CT/i/i = crMM(2/)/*o(2/ ~ yoikre)), 1^' = x,y. Using Eq. for Vg > s we obtain 

axx-ayy- • (i-^j 

Eq. (|l3|) coincides with jdW/Ey, given by Eq. (32) of Ref. ||l^a, for Ey 0. For Vg < s the contribution to iixx is 
exponentially suppressed and has an activated behavior |l^. For Vg ^ we have axx = e^foc'fc^T^/\/27r^/^?i^Wg, if 

1 > kBTlo/Uvg > s/V2vg, and Eq. (|l3|) if 1 > kBTlo/Uvg < s/V2vg . 

Notice that the terms in Eq. ( [lO| ) related to <Jyy{y) are totally absent in Eq. (|l2|). Now, in Eq. ( p^ we have 
\y~y'\ ~ ^0 and \kx\to ~ 10~^; then we can make the approximation A'o(|fc^||y — y'|) ~ ln{2/\kx(.Q\) — ^ — ln(\y — y'\/ to) 
where 7 is the Euler constant. The value of the integral overy'inEq. (^ is hi(2/|fc^|4)-7- (1/0F) ^^^dte-*^ ln\t- 
(y — yre)/^o|- For the gated sample and that with air above z = d, 4c?^ 3> ^q, the corresponding approximations in the 
long wavelength hmit are Rg w ln{2d/£o) — ln{\y — y'|/^o) and, with e 3> 1, -Ra ~ ln{2/kldio) — 2j — ln{\y — y'\/io), 
respectively. 

Now for r ^ we have axx 0, ayy — > and Eq. ( |lO| ) shows that /3(a), kx,y) behaves essentially as ^'□(y — yre)- 
Then the value related to the integral over t can be evaluated and gives a small contribution compared to that of the 
term ln{l/\kx\£o)- The result is 

1 3 f°° 

{-{cj ~kxVg) + S [In— — + -]} / dy p{uj,kx,y) = 0. (14) 

\Kxt0\ 4 

From Eq. (|lj) it follows that the EMP dispersion relation, with strong dissipation at the edges and for fc^T ^ hvg/£o, 
is given by {uj{kx) = uj) 

2 13 

UJ = kxVg + -[kxa°^ - ikldxx] [^'^j^-^ + 4]' (15) 

For V — 2 the EMP dispersion relation will again be given by Eq. ( p^ with the conductivity components multiplied 
by a factor of 2. In addition, because Vg has the ~ 1/B dependence, cf. Eq. (pi]), it will be multiplied, for v — 2,hy 
a factor of 2 if the edge field E^ is the same. As a result, ioY v = 2 the frequency w will be approximately twice larger 
than for v — 1. More exactly, the value of the ratio of these frequencies is 

uj{v = 2) _ ln2 

oj{v = 1) " 2Ml/|Mo|) + 3/4]^' 

where corresponds to = 1. For kxi^ ^ 10~®, the second term inside the curly brackets represents a 2% correction. 

Eq. ( p^ is valid for an ungated sample. If the sample is gated, repeating the procedure leads again to Eq. (|l|) with 
the factor [. . .] replaced by [^n(2d/€o|) + 2/7r]. If the gate is replaced by air, this factor is replaced by \ln(\ lk^dlQ\)\. 
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IV. NOVEL EDGE WAVES AT VERY LOW TEMPERATURES 



From Eq. ( |lo|) it follows that, even for T — s- 0, k^, y) — (w, kx)'^'^{y), y — y — 2/re: is only an approximate 
solution of this equation. A more accurate solution is obtained by the expansion 

OG 

p{iu, k,, y) = ^ k,)Hjy/£o) 

oo 

= ^V2^p("'(w,fc.)*«(y)*o(y), (17) 

where Hn{x) are the Hermite polynomials. Due to their orthonormality Eq. ( p7| ) is the exact expression for p(a;, kx,y)- 
Notice that this expansion is specific to the case when only the lowest LL is occupied. In addition, the terms 
n = 0,rt = 1,71 = 2, etc. correspond to the monopole, dipole, quadrupole, etc. expansions of p{Lu,kx,y) relative to 
y = yre- 

We now multiply Eq. (nn) by i?m(?//^o) a-nd integrate over y from y^e — Ay to yre + Ay. Taking into account that 
for very low temperatures^ hvg ^ ioksT) day^{y)/dy, (Jyy{y), and (Jxxiy) behave as '^^{y), we obtain 

°° 2"r?l 

^{u- kxVg)p^'^\u, kx) + (5 + m^) Y^{-—- f/^a„,,,{kx)p^^\uJ, kx) = 0, (18) 

n— 

where 

da: *,„(x)«'o(x) / da:' /Todfe^l k - x'|) *„(2;')^'o(a;') (19) 

oo J —oo 

and 5' = —4iayy/eiQ. Notice that for m — Eq. ( [l^ ) is equivalent to Eq. (|l2|) and correspondingly the terms 
related to 5-yy are absent. From Eqs. ( p7| ) -(|l9|) it follows that there exist independent wave modes, spati ally 
symmetric, p''{u!,kx,y), and antisymmetric, p'^'^Juj, kx,y), with respect to y = y^e- They are given by Eqs. ( |l7| ) 
and ( p^ with n even and odd, respectively. Notice that in Eq. (|l^) due to the assumption k^io <C 1 we can write 
Ko{\x\)w ln{2/\x\) — 7 for the ungated sample. For the gated sample we simply replace the kernel Kq in Eq. ( p^ by 
Rg « ln{2d/io) - ln{\y - y\/£o) and for that with air by Ra « ln{2/kld£a) - 27 - ln{\y - y'|/4)- 



A. Symmetric modes 



Considering only the term n = in the expression for p*(w, kx,y), Eq. (18) for m = gives 

[~{uj-kxVg) + SaQo\p^''\u,kx)=Q, (20) 



where amikx ) = -ln{\kx\lo) + 3/4. With this value of awikx) and p(o)(w,fca;) 7^ 0, Eq. (|20l) gives the dispersion 
relation (|5|). Because p^{oJ, kx, y) behaves spatially as ^^{y) in this approximation, we will refer to it as the dispersion 
relation of the monopole EMP. For the sample with a gate we simply have to replace ooo by Oqq = [/n(2d/^o) + 2/7r] 
in Eq. ( po| ) and for that with air by Ogg = ln{l / k'^i^d) . 

Corrections to Eq. (|2^) and further symmetric branches are obtained by keeping only the terms n = and n = 2 
in the expression for p^iuj, kx, y) which gives 

p^Lo, kx,y) = p("Hc^, kx)^l[y) + 2x/2p(2)(c^, kx)-^ o{y) ■ (21) 
From Eq. ( p^ ) for m = we obtain 

[-{uj - kxVg) + S aoo]p^°^(w, kx) + 2V2S ao2P^''\uj, kx) = 0, (22) 

and for m = 2 

[-{lo - kxVg) + {S + 2S') a22]p(') (^, fc.) + [{S + 2^')/2\/2]ao2P^°^ (w, fc.) = 0, (23) 



where we write amn{kx) = o^mn in order to simplify the notation. For a nontrivial solution of the system of Eqs. ( |22 
and ( p3| ) the 2x2 determinant of the coefhcients must vanish. This gives two branches u)'!^{kx) and ui^_{kx)- For 
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\kx\(o ^ 1 a numerical evaluation gives 002 = —0.353, 022 = 0.250, and = 1/8. All values remain the same 
for gated samples or those with air except aoo which changes as indicated above. If we neglect the coupling terms, 
by formally setting 092 (fcx) = 0, Eq. (|2|) gives the monopole EMP dispersion relation ( [l5| ) and Eq. ( p^ ) the pure 
quadrupole EMP dispersion relation 

uj = k^vg + (5 + 25')/4 

= k^vg + ^ [k^al^ - ikld.,^ - ■ (24) 

If we neglect k^Vg and the dissipative terms, Eq. ( |2^) takes the form of Eq. (14) of Ref. |ll] for the j = 4 branch 
which has 5 charge oscillations. As it stands Eq. (|24|) corresponds to only 3 oscillations. The difference is to be 
ascribed to the very different density profile used in Ref. [lll| for a compressible liquid in a very wide strip. Notice 
that Eq. ( p4| ) is valid for samples with gate or air as well. 

For finite ao2 the two branches resulting from Eqs. (|2^ ) and (^3|) are given by 

= kxVg + (l/2)[S'(aoo + a22) + 25'a22] 



^± — r^x^Jg 

±(l/2)^[5(aoo - 022) - 2S'a22? + 45(5 + 2S')al^. (25) 

If we put ao2 = 0, i.e., if we neglect the coupling between the branches, then the uj'L{kx) branch is given by Eq. ( [2^ ) 
and the u;^(fc2,) branch coincides with Eq. (p^). It can be shown that the term oc 0^2 under the square root sign is 
much smaller than the other term. Then from Eq. ( ^5|) we obtain 

S{S + 2S')al 
SaoQ — {S + 2S')a22 ' 

and 



0.; = k^vg + Saoo + . ^ 7. ^ (26) 



s 7 I (-c I oc'^ S{S + 2S')aQ2 

w_ = kxVg + [S + 2S )a22 - , , . (27) 

baoo - (o + 2b')a22 

Further, for these very low temperatures we can distinguish between (i) strong dissipation, for which k^cTy^ <C 
Adyy/ig aucl (u) weak dissipation, for which kxUy^ 2> Aayy/lQ. Notice that the damping of the purely quadrupole 
EMP-Eq. (E4)- is such that in case (i) we have ^ \^u>\ whereas in case (i) the opposite inequality holds. The 
damping of Eq. ( p^ ) is determined by the dissipative conductivities <Jyy{y) and (Jxx{y)- The two contributions differ 
by a very small factor k'^l'^. As a result, the damping of the wave, oc k'^axx, can be usually neglected. Notice that 
Eqs. (p^)-(p7|) are valid for gated or "air" samples as well with ooo replaced by or agp. 

For definiteness in numerical estimates we will use parameters pertinent to GaAs/AsAlGaAs heterostructures. 
As will be shown below, both cases (i) and (ii) are experimentally realized depending on the values of Vg and T. 
For Vg ^ s the damping in Eq. ( p4[ ) is exponentially suppressed. The condition Vg < s requires a smooth energy 
dispersion near the edges. This possibility exists in the Hartree approximation for the confining potential but not 
in the Hartree-Fock approximation where the exchange leads to a logarithmically divergent Vg ]1^. However, when 
correlations are taken into account, a smooth energy dispersion results near the edges and Vg is small ^M. In GaAs- 
based heterostructures the most common case is Vg > s = 2.5 x lO^cm/sec. In this case using Eq. (jlj) for 1^=1, 
c' = h{ehi4)'^ /2pvs, hi4 = 1.2 x 10^ V/cm, and pv — 5.31 gm/cm'^, we obtain 



P = O.ief 3^-3/2^-4^ (28) 



ioa^yx 

where T — T/1°K,B = -8/(10 Tesla), and Vg — Vg/s. Eq. ( p8| ) is valid ioi v — 2 as well because of the scaled 
quantities. For f = 1,5 = 1, and ig = 2 Eq. (|2|) gives CTyj;/4cr°^ = 10"^ ^he estimated field leads, 
for _B = 1, to « 4 X lO^cm/sec > s. Further, if we assume that ^yyjlQ gives approximately the value of the 
diagonal conductivity in the edge strip of width 1^ « 8OA, then, because for strong magnetic fields (lOcT* ^ 1) 
UyyllQU^^ ~ XjijJcT* = 10~^, we obtain an effective scattering rate 1/r* « ujcdyy/ioay^ « 2.6 x 10^^/sec in this strip. 
This is approximately ten times larger than the scattering rate for a mobility p — 10^ cm^/V sec. 
From Eq. (|2^ ) we obtain 

P^'^^l_,f-'B-"V,lkM. (29) 

1 3a; I ^(Tyy/lo ^ 
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where we assumed again that the term Vgk^ can be neglected. Then, for T = 1, B = 1, and Vg = 2, the RHS of Eq. 
(^) is approximately equal to 25|fca;|£o- Only for 1 ^ \kx\io > 4 x 10~^ does the quadrupole EMP become weakly 
dissipative. For B — 1, lower temperature f — 0.1, and steeper confinement Vg — 4, Eq. ( p9|) gives 1 3? w | / 1 3[j | « 
5 X 10^|fca;Ko- In this case the quadrupole wave is very weakly damped for 1.25 x lO^cm"^ \kx\ > 2.5 cm~^. In this 
region the implicit low-frequency condition \uj\ <C oJc is well satisfied since it corresponds to \kx\ <SC 1.6 x lO^cm"^. 



B. Edge helicons 

We now analyze further the general formulas of subsection A. We first assume that cTy^lkxlK ^ ayy/t^. Then Eq. 
gives 

ut = kxVg + \{S + 2S'){l-^), (30) 

where K — oqo — 1/4 = ln{l/\kx\to) + 1/2; for gated or "air" samples ooo in K is replaced by Oqq and agp, respectively. 
Because ln{l/\kx\lo) ^ Ij we see that the coupling with the monopole EMP does not change the dispersion almost at 
all as compared with that given by Eq. (HJ). As a result the quasi-quadrupole EMP, described by Eq. (|30|), is weakly 
damped for cr°^|fc2;| > AcTyy/il and strongly damped for AcFyy/lQ > (Jyx\kx\ ^ ^yy/i^K. We call the wave described 
by Eq. ( |30| ) quasi-quadrupole EMP because it follows from Eq. (|2^) that 

Qj^, = -2V2ao2i^ = if » 1; (31) 

that is, p^^\ui, kx) is the dominant term on the RHS of Eq. (^T|). The same holds for the "air" sample. However, for the 
gated sample the monopole and quadrupole terms are comparable if Kg < 3. The condition of very weak damping for 
the wave (|3(]| ) can also be expressed as wir* 3> vrc/iT where tq = e^/ehcjcio- For v — 1 we typically have vq 1. This 
condition resembles the high-frequency limit used in Ref. [Q. However, here r* is related to dissipation processes only 
near the edge. In addition, in contrast with Ref. ||l[ we consider an essential decrease of the conductivity components 
and of the electron density over a finite length Ay from the edge. For very low temperatures, ksT <C hvg/io, we have 
Ay « ig which is much smaller than the length over which the density nQ{y) decreases substantially in the model of 
Ref. P . 

FoT^xlkxlK > ayy/el Eq. (P gives 

As can be seen, taking into account the coupling with the quadrupole EMP changes the phase velocity of the monopole 
EMP by a very small amount (< 0.1%) but it gives a principally new contribution to the damping in comparison 
with Eq. (p^. Now the lu^ branch, for the typical EMP condition \kx\io ~ 10^^, has a damping c>c ayy/i^K which 
is much stronger than that oc axxk'^K of the pure monopole EMP. 

The wave described by Eq. ( ^ ) can be called a quasi-monopole EMP because it follows from Eq. (|2^) that 

P^°\^.kx) ^ 16V2ao2a°xk.K 

pi^){u:,kx) alxkx-^^Oyy/^l' ^ > 

and, due to K 1, we have |p^°)(w, fca;)/p^^^(w, fc^;)] '> 1. Now, for weak dissipation {cfyxkx ^ ^^yy/^o) have 
p(o)(w,fc^)/p(2)(t^^fc^) ~ _8if and for strong dissipation (cr°^fc^if > ^yy/el > cr^^kx/A) p^^\uj,kx) / p^'^\uj,kx) « 
—2 {{ay^kx^"^/ cTyy)K . Thus, if the phases of the two components are shifted by tt in the first case, in which we call 
the branch described by Eq. (P2| ) a modified monopole EMP (MMEMP), they are shifted by 7r/2 in the last one. 
This last case corresponds to the frequency region w^t* ^ vr^jT^ ^ {ujI^t* ji^K + 1)) and the frequency can 
still be considered as high compared to 1/r* . In this frequency region we call the cj^ branch the high-frequency edge 
helicon (HFEH) and denote it by In this region, due to the almost 7r/2 shift between p^^^ and we obtain the 

following remar kable property of HFEH described by Eq. (||). If the HFEH charge oc Re{p'{ujl^,kx,y)R+{x,t)\ 
along y has a pure quadrupole character cx \p^'^^ \ for some phase of the running wave i?+(x, t) — exp[—i{uj^f^t ~ kxx)], 
then after approximately a ±7r/2 shift it acquires a pure monopole character oc Ip^*^-*]; p^{cj^§, kx,y) is given by Eq. 
(pT|). This HFEH shows three charge oscillations along y whereas the relevant branches of Refs. j^] and |jll|, with 
Reoj oc kxln{l / kx) resembling 3?ijj|(^ of the HFEH, show only one oscillation. 
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Notice that 3?w^, given by Eq. (||) is independent of T whereas 3u;^ cx or if 1 > kBTio/hvg > s/V^Vg. 
That is, in contrast with Ref. [0 we find that the damping of the MMEMP and that of the HFEH scale with 
temperature and are not quantized in the QHE plateaus. In addition, these waves have a characteristic length £q 
which is different than the length £y of Ref. Q; also, the term oc crjj^ is different than the that of Ref. in the factor 

containing the logarithm. In addition, £q{i' = 2)/£o{i> = 1) = \/2 here whereas £v{i' = 2)/£v{v = 1) = 4 in the high 
frequency limit of Ref. Moreover, for e ^ oo it follows that £-u{v — ^) ^ Q whereas £q is independent of e. 
We now consider the case of very strong dissipation 

~^yy/^l » ^IMK. (34) 

Then Eq. ^ gives 

u'_ = k,Vg + L [ial^k, ~ '^l^yy/f-l] ■ (35) 

This is again a quasi-quadrupole wave since |p'^'(w, k^) / p^^'' {oj , fcx)| ~ ^yy I '^yx\kx\£o ^ 1- Also, although ^uj^/kx is 
essentially larger than in Eq. (|30|), Eq. ( |3^ ) gives an aperiodic damping, j^ftwlj <C ISJo-iil. Further, assuming that Eq. 
( ^ ) is valid, we obtain from Eq. ( p6| ) the dispersion relation of a new wave that we call low-frequency edge helicon 
(LFEH) as 

= kxvg + S{K - i) - -^^(i^ - \). (36) 

Despite the strong dissipation condition (|3^), which entails w^l^r* ^ frQ/n < 1, the LFEH is very weakly damped 
since |3?a;|^|^| ^ jO'ix'^IJ^I. The frequency range of the LFEH is similar to the low- frequency limit of Ref. M but here 
T* is related to strong dissipation processes only near the edges. Also, ^uj^^/kx differs little from ^uj^^/kx or that 
of Eq. (p^ . However, the damping of the LFEH oc {ayxkx£o)'^ln{\ / kx£o) / cFyy has a very different form than Sw^jj' 
of Eq. @ or 3cj of Eq. (|l|). 

We further notice that, in contrast with Ref. [Q, the real part '^ujj^^ is independent of temperature whereas the 
imaginary part i. e., the damping, is not quantized and varies as or T^^; the behavior occurs if 

^ and 1 ^ kBT£o/hvg ^ s/V^Vg. That is, the LFEH has a characteristic length very different than £y in the 
low-frequency limit. It follows from Eq. (^) that 

P^°H^j^^,,k^ 



-4V2ao2 = 2. (37) 



With this result and Eq. ([2l| ) we obtain the dimensionless charge density profile of the LFEH, pEH{y) — 
^£op' (uj^'h ,k,,y)/p^°^ {LJ^fi , fcx ) , as 

PEHiy) = V^^o[*o(^) + V2^2{y)My)]- (38) 

In Fig. 2 we show PEH{y) (solid curve ), its monopole component (term cx ^q, short-dashed curve), and its 
quadrupole component (term (x ^9^2, long-dashed curve ). For contrast the dotted curve represents the normalized 
unperturbed electron density no(?/)/7io- As can be seen the monopole and quadrupole contributions arc of the same 
order of magnitude and the resultant pehIv) has an oscillatory behavior with two oscillations, one to the right and 
one to the left of the edge at y = y^e. This is in sharp contrast with the "usual" EMPs of Refs. [Q and the j = 
mode of Ref. 

C. Antisymmetric modes 

Considering only the term n = 1 in the expression for p°''^{uj, fc^,, y) and using Eq. ( p^ for m = 1 we obtain 

[~{u~h,Vg) + {S + S')aii]p^^\u,kx)^Q- (39) 

the numerically obtained value of an is 0.5. For p(i)(w,ft^) ^ Eq. (||) gives the dispersion relation of the pure 
dipole EMP as 

uj = k^vg + (l/e)[fc.a° - 2idyy/£l]. (40) 
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If we neglect kxVg and the dissipative term, En. (^) will take the form of Eq. (14) of Ref. [|T| for the j = 2 which 
shows three charge oscillations whereas Eq. (40) corresponds to only two oscillations. However, in contrast with Ref. 
[ pl| , besides the term k^Vg and the microscopically treated dissipative term, the Hall conductivity cr"^ is quantized 
for V = 1 01 2. 

Corrections to Eqs. (^) and (^) and further antisymmetric branches are obtained by keeping only the n = 1 and 
n = 3 terms in the expression for p°'''{uj, kx,y). Then 

p'^^(L^,fc,,y) = V2p(i)(c^,fc,)*i(y)«'o(y) + 4V3p(3)(^,fc^)5-3(y)^-o(y). (41) 
From Eqs. (|l^ ) and (|9|), for m = 1, we obtain 

[-{lu - k^vg) + [S + S') an]p''^HL^,k^) + 2V6iS + S') aisp^^^ (w, fc.) = 0, (42) 

and for m = 3 

[-{uj - k,vg) + [S + 35') a33]p^''(^, + (1/2V6)(5 + 35')«i3P^'^(^, K) = 0. (43) 

Again the vanishing of the 2x2 determinant of the coefficients gives the two branches cj^i" and u}";^ . In the long 
wavelength limit we numerically evaluate ais = —0.204 and 033 — 0.166 =1/6. If we neglect the coupling between 
the modes, by formally setting ai3 — 0, Eq. (42) gives the dispersion relation for the pure dipole EMP, Eq. (^), and 
Eq. (H) that for the octupole EMP 

0J^k,Vg + (l/3e) [fc,a°, - Qidyy/ll] . (44) 

If we neglect k^^Vg and the dissipative term Eq. (^}takes the form of Eq. (14) of Ref. jll] for the j = 6 branch that 
shows seven charge oscillations. As it stands, Eq. (ij) corresponds only to four charge oscillations. 
For finite 013 the two branches resulting from Eqs. ( p^ ) and (^3|) are given by 

c^l^ = k^vg + (l/2)[5(aii + 033) + 5'(aii + 3033)] 

±(l/2)y^[5(aii - 033) + S'{aii - 3a33)]2 + 40^3(5 + S'){S + 35'). (45) 

If we set ai3 — 0, uj'^ and u}°f give the dipole and octupole branches given above by Eqs. (^) and (^), respectively. 
For weak dissipation we have kxCfy^ ^ '^^yy/^o neglect damping, we obtain 



and 



kxVg + i6/5e)a°^kx (46) 



kxVg + {l/7e)a}'^kx. (47) 



We call the waves corresponding to Eqs. (g§) and (g^) modified dipole (MDEMP) and octupole (MOEMP) EMP, 
respectively. As we now show, some of their properties are essentially different than those of the pure dipole (Eq. 
(|g)) and octupole (Eq. (||)) EMPs. 

Then from Eq. (|2|) for the MDEMP we calculate 

4JJ^^^=20V6ai3 = -10.0. (48) 

With this ratio and Eq. 

(|l|) the dimensionless charge density profile p"^ {lu%' , k^ , y) = ^/niop'"'{LJ± , kx,y)/ p'^^'>{uj± , kx,y) of the MDEMP 
takes the form 



PMDiy) = pXiy) = V27r£o [*i(y)*o(y) - (l/2)*3(y)f o(y)]. (49) 

In Fig. 3 we show PMoiy) (solid curve), its dipole component (term cx ^'i(y)5'o(y), short-dashed curve), and its 
octupole component (term cx '^3(y)^>Q{y) long-dashed curve). Again the dotted curve represents the normalized 
unperturbed electron density no(y)/no. As can be seen the dipole and octupole contributions to pMoiy) are of the 
same order of magnitude though uj'^^ = lomd-, given by Eq. (|46|), is a bit different (< 20%) than SRw of the pure dipole 
EMP given by Eq. (Eoh. As seen in Fig. 3, the MDEMP has four charge oscillations whereas the dipole EMP has 
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two. Thus, the corresponding density profiles are qualitatively different although the phase velocities are close to each 
other. 

For the MOEMP, given by Eq. @, we calculate from Eq. (|3|) 

4^^^^^ - 2.109/2V6 « 1/%^; (50) 

the corresponding result for pMoiv) = P°^{y) is 

PMoiv) = \/2^4 [^limoiy) + 2«'3(y)*o(y)]. (51) 

In Fig. 4 we plot the same quantities as in Fig. 3 but now for pMoiv)- As can be seen, the spatial behavior of pMoiv) 
is quantitatively different only than that of the pure octupolc EMP. The phase velocities of these two EMPs, as follow 
from the dispersion relations, are substantially different. Notice that for strong dissipation, kxCfy^ ^ 4(Tj,j,/£q, both 
Lo'^ branches are strongly damped. 



V. DISCUSSION AND CONCLUDING REMARKS 



We have introduced a realistic model for the confining potential Vc{y) and treated mainly the case where v = \ m 
the interior part of the channel and the formation of dipolar strips [|l8| near the edges is impossible in the assumed 
QHE regime. We have taken Vc{y) sufficiently steep at the edge that LL flattening can be neglected [Q- |2^. As 
for V — 2, we have neglected the spin splitting. This is a reasonable approximation in the bulk of the channel but its 
validity near the edges is not clear in view of the work of Rcfs. jl^] and j2l]. Though we have used a simple form for 
Vc{y), near the edge the results remain valid for potentials Vc{y) of different form that are smooth on the scale of ^o- 
In such a case, e.g., the last term on the RHS of Eq. (||) should be replaced by Vc{yoa) and Vg will be given again 
by Eq. ( |ll| ) with kre determined by foikre) = 1/2. Thus, Vg is the only parameter related to the form of Vc that 
influences these edge waves. 

For a comparison of the theory with the experiment the details of the sample geometry are necessary but often 
they are not given. For instance, the thickness dg of the sample is not given in Ref. it is also not given neither in 
Ref. , where the bottom of the sample is metallized, nor in Ref. Now concerning Ref. |^ it can be seen that 
ds = 400 pm for the sample with perimeter P = 1.7 cm and the observations were made at T = 0.4K and B — 13 
Tesla {u = 1). As we now show, the first three edge helicon modes k^o = 27r/P, k^i = 2kxQ, and kx2 = Sfc^jo do not 
lead to an equidistant spectrum. Indeed, in addition to the case of the sample with air treated above, there is air for 
z < —dg (a foam plastic base with small dielectric constant). Then for 2kxds <C 1 we arrive at Eq. (^) with e « 1, cf. 
Ref. |2^, which should be substituted in all formulas involving a homogeneous e; for 2kxds 3> 1 we obtain the case of 
the "air" sample treated above. Taking into account that 2kxods = Airdg/P ^ 0.3, 2kxids ~ 0.6, and 2kx2ds ~ 0.9, 
we can explain qualitatively the observed, not equidistant spectrum. A quantitative comparison is impossible because 
the condition |A;j;|iy 1 is not satisfied (|/ca;|M^ 1). 

The sample of Ref. |^ is a circular mesa with diameter D = 540 pm, height ~ l^m, in the middle of a wide 
GaAs/AlGaAs chip with thickness ds — 500 pm. The condition 2|fc:i,|rfs ^ 1 is well satisfied even for kxo = 2/D = 37 
cm~^; the minimum value of kx, kxmim involved in the experiment satisfies kxmin 3> kxo- As a result, the conditions 
"^kxminds ^ 1 and kxminW ^ 1 are well satisfied as well. Because the square "pulser" gate, of width Lp — 10/im was 
much smaller than the circumference ttD, the initial charge distribution can be assumed to have a rectangular form 
and therefore to have an essential contribution from the kxn = ^nkxo, n = 1, 2, .., modes distributed in the interval 
kxO ^ \kxn\ < T^/Lp. It is natural to assume that the typical kxt — \kxnt\ ~ '^/'^Lp] for a rectangular form of the 
charge distortion along x, we have approximately a 50% contribution to the total spectral density for \kxn\ < kxt- 
The model of the "air" sample, described briefly in the Appendix, fits perfectly to the experiment ||]. The distance 
d at which the 2DEG is situated beneath the surface is not given but from the mesa height it can be estimated to 
be d ~ lOOOA. This gives Ka{kxt) ~ 10.5 for B ~ 5.1 Tesla. Further, assuming Vg = 2, i.e., = 5 x 10^ cm/sec, 
using B ~ 0.51. T — 0.3, and v ~ 1, in the middle of a wide chip, obtained from the parameters of the experiment, 
and using Eq. (|2^), we obtain <7yy/£oay^ « 7.5 x 10^''. It follows that kxdoKa{kxt) ~ 1.6 x 10~^ ^yy/^o^^yx which 
corresponds to the first condition treated in Sec. IV. B. 

Notice that here kxt^^o ~ 1.5 x 10"'^ < '^o-yy/ i^ay^ w 3 x 10"'^, i.e., we are dealing with the case of strong dissipation 
as detailed after Eq. (|27|) . It follows that all modes, apart from the HFEH, are strongly damped. For instance, both 
the pure dipole (Eq. po|)) and quasi- quadrupole (Eq. (pO|)) modes have Slj « 3.4 x 10^°/sec which is a very large 
value compared to the damping rate of the experiment '^uj « 2 x 10^/sec for v — 1. Moreover, during the observed 
period for the travel around the edge Tg « 3.5 x 10^^ sec the amplitude of the quasi- quadrupole EMP or dipole EMP 
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should practically vanish due to the exponentially small factor < exp{—100). We are left with the high-frequency 
edge helicons described by Eqs. ( |32|) and (|3^). Eq. (|3^ ) gives a decay rate Qujeh ~ 1-9 x 10*/sec which is in 
good agreement with the observations |^ . The corresponding group velocity for the HFEH, obtained from Eq. (|3^) , 
is Vg{kxt) ~ (2/e)<7y2;[^"-(l/^L^o'^) — 2] and gives a period T = ttD / Vg{kxt) ~ 3.4 x 10~^/sec which is in excellent 
agreement with the experimental value. Because after the first trip the pulse became 700 ps wider, we can estimate 
the range of lS.kx, around kxt, which gives the most essential contribution to the pulse. Subtracting Vg{kxt + 
from Vg{kxt — ^kx/2) leads, after some calculations, to /S.kx/2kxt ~ 0.4. This then gives approximately kxmax — ^■4:kxt 

and kxmin — O.G/Ca;^. 

Another important ingredient of our theory are the calculated damping rates. As shown above, they agree well 
with the observed rates The effective r* for the work of Ref. is approximately 2.5 x 10"^"'^ sec. This is several 
decades shorter than that extracted from QHE measurements which is of the order of 10"'^ sec |Q. The difference is 
to be ascribed to the fact that in our model the dissipation is localized near the edges whereas in that of Ref. |^ the 
dissipation occurs throughout the channel homogeneously. The latter is a reasonable assumption for relatively high 
temperatures which have not been considered in the present work. 

Finally, it is worth noticing that in our microscopically calculated dispersion relations the quantized Hall conduc- 
tivity ayx appears explicitly. Though not shown here graphically, this accounts for the existence of plateaus in the 
transit times of the signals, as a function of magnetic field, observed recently in Ref. |^ and accounted for by the 
replacement of the quantity by ayx- The present theory holds when only n=0 LL is occupied. The coupling between 
edge excitations of different LLs will certainly affect the EMP modes presented here and an appropriate extension of 
the theory is being planned. 
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APPENDIX 

As a model of a gated sample, we consider a 2DEG at z = with a metallic gate placed at a distance z = d away 
from it and with a dielectric constant e for z < d. Taking the Fourier transform, for z < d, with respect to x,y, and 
t, of Poisson's equation for the time-dependent charge density p{x, y, z, t) gives 

^2 

e[^^ - kx, ky, z) = 4:TTp{uj, kx, ky)S{z), (52) 

where fc^ = fc^ + ky . For z < we have 

(t)-{uj, kx,ky,z) = A{uj, kx, ky) e*"^, (53) 

and for < z < d 

(/)+(w, kx, ky, z) = B{uj, kx, ky)e'''' + C(w, kx, ky) e'''\ (54) 

Two boundary conditions are (j)j^(ui,kx,ky,d) = and <p^{uj,kx,ky,+0) = <p^{uj,kx,ky,^0). Integrating Eq. (|5 
from z = —0 to z = +0 gives the third condition 



b+{uj,kx,ky,z) d(j3-{iv,kx,ky,z) a f u i, \ 

el U=+o - \z=-o\ = -^T^P(^,kx,ky). (55) 



From Eqs. (|53|)-(|55|) and the first two conditions we obtain 

cj}{ij, kx,ky, z = 0) - 2^p{uj, kx, ky){l - e-2fe'i)/eA: (56) 

which gives 

(j}{u;, kx,y) = — / e'''yy(t){uj, kx, ky, z = 0)dky 

2 

e 



[Ko{\kx\\y~y'\)~Koi\kx\V¥^7F+^M^,kx,y') dy . (57) 

Notice that for d ^ oo Eq. ( |57| ) coincides with Eq. (||). 

If the gate is replaced by air, for z > d, a similar calculation leads again to Eq. (^) with the kernel replaced by 
Ra - Koi\kx\\y - y'\) + [(e - l)/(e + l)]Koi\kxWiy - y')^ + M^). 
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FIG. 1. Unperturbed electron density no{y), normalized to the bulk value no, as a function of y/£o measured from the right 
edge taken as the origin. The solid and long-dashed curves are obtained from the models of Refs. jl| and [Q, respectively, as 
explained in the text. The short-dashed curve is the profile of the present work. 

FIG. 2. Dimensionless charge density profile pEH{y) of the low-frequency edge helicon (LFEH) described by Eq. (^) as 
a function of y/lo (solid curve). The dotted curve represents no{y)/no. The dashed and long-dashed curves are the pure 
monopole and quadrupole contributions, respectively. The oscillatory behavior of PEniy) near the edge is in sharp contrast 
with the "usual" EMP of Refs. [§ and [|ll[. 

FIG. 3. Dimensionless charge density profile pmoiy) of the MDEMP described by Eq. ( p9| ) as a function of y/lo (solid 
curve). The dashed and long-dashed curves are the pure dipole and octupole contributions, respectively. The dotted curve 
represents no(y)/no. 

FIG. 4. As in Fig. 3 but for pMo{y) as described by Eq. (|l|). 



12 



no 



1.25 



1.00 



0.75 



0.50 



0.25 



0.0 




-20.0 -15.0 -10.0 -5.0 



y/i 







1 




2 




3 




4 



